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^ ■ Abstract 

VO . The SNla data, although inconclusive, when combined with other observations 

makes a strong case that our universe is presently dominated by dark energy. We 

f^ ' investigate the possibility that large distance modifications of the curvature of the 

universe would perhaps offer an alternative explanation of the observation. Our cal- 
culations indicate that a universe made up of no dark energy but instead, with a 

i-^ \ modified curvature at large scales, is not scale-invariant, therefore quite likely it is 

ruled out by the CMB observations. The sensitivity of the CMB spectrum is checked 
2 . for the whole range of mode modifications of large or short distance physics. The 

^ \ spectrum is robust against modifications of short-distance physics and the UV cutoff 

Cd ■ when: the initial state is the adiabatic vacuum, and the inflationary background space 

is de Sitter. 
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1 Introduction 

Based on the theoretical cosmological models of inflation, the interpretation of the current 
astrophysical observations such as the SNla data [|l|, suggest that our universe contains a 
large amount of dark energy |Q . 

However, alternative models, free of dark energy, which may fit in the allowed range of 
parameters suggested by observation, are not excluded. In this paper we investigate claims 
to a possibly different interpretation of the SNIA data, for these alternative cosmological 
models: a FRW universe with no dark energy but with a modified curvature at large 
enough distances. The hope then is that either the Friedman equation for the expansion is 
modified, or that the light from SNla that reaches us, while passing through these regions 
of different curvature, would be deflected, thereby "appearing" to have the same effect as 
an accelerating universe[}|. 

We examine metric perturbations in this modified background geometry (traced back at 
the time of inflation). Metric perturbations are responsible for the generation of the large 
scale structure and temperature anisotropics of the CMB. The inflaton field (in 4 dimen- 
sions), through the Friedmann equation, determines the expansion rate H for the curvature 
of the background geometry. The metric perturbations satisfy a Klein-Gordon scalar field 
equation, minimally coupled to gravity [^ . The scalar field has a generalized mass squared 
Vtn{TiY that receives the contributions of two terms: the field frequency squared and the 
field coupling to the background curvature term. The coupling of the field to the curvature 
results in a modified propagation at long wavelengths since the curvature of the universe 
is modified at large distances compared to the intermediate scales. Examples of modified 
gravity can be found in [|, H, |^, P] . Then, the modified propagation of wavelengths of the 
same scale as the background curvature deviation scale can be attributed to a nonlinear 
dispersed frequency of the field at those wavelengths, for as long as the generalized mass 
squared, VLn{riY in the field equation, remains the same. This equivalence noticed in is 
very useful for calculating the effects of modified large distance curvature in observations. 

Our model consists of a (2-parameter) family of nonlinear dispersion relations for the 
generalized frequency of the field, that take account of the modification of the curvature 
at large distances. The family of dispersion functions is nearly linear for most of the range 
k < Mp, except a nonlinear deviation centered around some low value of momenta ko. It is 
this deviation bump that reflects the modifications of the generalized frequency of the field 
at low momenta k^ due to the modification of the curvature at large distances k^^. The 
dispersion function introduced in Sect. 2, although nonlinear in the transplanckian regime, 
it is nevertheless a smooth function there, asymptotically approaching a constant value at 
time- infinity, thus having a well defined initial vacuum state |T^. The analytical calculation 



of the CMB spectrum is based on the Bogoliubov coefficient method. The details of the 



exact solutions for this class of dispersion functions \TU\ are given in the Appendix. 

In Sect. 3 we check the sensitivity of CMB spectrum to the bump parameters ko and B 
(scale location and amplitude) that control the deviation behavior from a linear frequency 

^We thank A. Riotto for bringing this idea to our attention. 

^See however Ref. H for constraints on models with spatial variations of the vacuum energy density. 



dispersion at low values of the momenta; i.e., the allowed range of curvature modifications 
at very large or very short distances that may agree with observation. We use CMBFAST 
in Sect. 3 to plot the spectrum, by replacing the standard primordial power spectrum 
S%{k) with that derived analytically in Sect. 2 for the model considered. We comment 
and summarize the results in Section 4. It is shown that the CMB spectrum is sensitive 
only to the choice of the initial vacuum state and the departure from linearity in the low 
momenta regime. However, for an adiabatic initial vacuum state, the CMB spectrum of 
a de Sitter expansion does not depend in the details of nonlinearity in the transplanckian 
regime 
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2 The Model 

The generalized Friedmann-Lemaitre-Robertson- Walker (FLRW) line-element in the pres- 
ence of scalar and tensor perturbations, takes the form []TB| 

ds'^ = a^{ri) {j-drf + [5ij + h{rj, n)QSij 

+hi{r], n)^ + hg^ir], n)Qij dx'dx^ i , (1) 

where r] is the conformal time and a{ri) the scale factor. The dimensionless quantity n is 
the comoving wavevector, related to the physical vector k by k = n/a{ri). The function 
{h, hi) and hg^ represent the scalar and tensor perturbations respectively. 

The power spectrum of the perturbations can be computed once we solve the equations 
in the scalar and tensor sector. The equation for the metric perturbations corresponds 
to a Klein-Gordon equation of a minimally coupled scalar field, /i„, in a time dependent 
background^ 

/i'^ + n„(r7)Vn = 0, (2) 

where the prime denotes derivative with respect to conformal time 77, and the generalised 
comoving frequency isQ 

a,(r]Y = n^-- = a''e--. (3) 

a a 

The dynamics of the scale factor is determined by the evolution of the background in- 
flaton field </>, with potential V{(j)), and the Friedmann equation. There are mechanisms 
that may produce different scale factors by modifying gravity at large (e.g. §, |^, 0, §]) or 
short distances (|13[).The present large distance modification scales can be traced back in 
time and would correspond to deviations in the primordial scale factor and spectrum. We 
can denote this "distance-dependent" scale-factor by A. The coupling of the field to this 
background curvature results in a modified propagation of the field at long wavelengths. 
Therefore, modifications of the scale factor or curvature (A) of the universe at large scales 

^We refer the reader for the details of the procedure to Refs. U% and related references H|. 
*Note that from here on we use the symbol a instead of a{ri) for the scale factor. 



can be attributed to a dispersed effective frequency {nejf), sucli tliat the generalized co- 
moving frequency Eq. (Rl) remains tlie same, in tlie following manner 






n^'if = "' - ^ = "h, --. (4) 



rieff denotes the dispersed comoving frequency of the field due to absorbing the modifica- 
tion terms to the curvature, A" /A. Therefore, the dispersion function for the generalized 
frequency results from the modified curvature at very large and very short distances. It 
deviates from linearity at small momentum k and asymptotically approaches a constant 
value in the transplanckian regime. 

The dispersion relation for the generalized comoving frequency ^n{v) is simplyQ [p!0[] : 
^n{v) = a{ri)F[n/a{ri)]. The 2-parameter family of dispersion functions F[k] of our model 
(see Fig. [^) is: 

F[k]^ = (e-kl)Voix,Xo) + k^Viix-Xo) + kl, (5) 

^o(^'^o) = rT^+(i + e-)(l + e(— 0))' (6) 

Vlfx — Xo) = —B— -, r-TT , (7) 

where the dimensionless wavevector is x = k/kc, kc = Mp is the cutoff scale, feg <^ kc, (i.e. 
Xq ^ 1) is the value at which we deviate from linearity at low momentum, the deviation 
amplitude is controlled by the parameter B, and the constant parameter ki < kc is the 
asymptotic value of the frequency in the transplanckian regime {k — ^ oo). C,E,B,xq are 
dimensionless parameters. 

As already discussed in Refs. [jll], [l^, 0, Eq. (Q) represents particle production in a 



time-dependent background |T^, 19|. We will follow the method of Bogoliubov transfor- 
mation to calculate the spectrum. The frequency Qniv)"^ (which is the same as a 'time- 
dependent mass squared' term) goes asymptotically to constant values at late and early 
times. Therefore the initial and final vacuum states are well defined. At early times the 
wavef unction should behave as a plane wave: 

u -^ ^ ^-i^nV (Q) 



But at late times one has a squeezed state due to the curved background that mixes positive 
and negative frequencies. The evolution of the mode function yU„ at late times fixes the 
Bogoliubov coefficients a„ and /?„, 

^„ ^,^+00 -^e-^""*" + ^^e+^^""'' . (9) 

with the normalization condition: 



2fi°«* A/2fi°«* 



\an\^ - \Pn\^ = I . (10) 

^From here on, we absorb the term a" /a of Eq. (|j) into the definition of the dispersion function F[k]. 




Figure 1: Shown is one of the members of the family of dispersion relations il„(r/)'^ as a function of: 
(a)fi^ = a'^F[k]'^ as a function of the physical momentum x; (b) conformal time rj. The variable x = k/kc 
has been shifted by 1 such that the regime of linear dispersion relation corresponds to small positive values 
of momenta, i.e., x <S^ 1. 



In the above expressions, fi^" and fi""* denote the asymptotic values of r2n(^) when rj — > 

=Foo. 

Details of the exact solution for Eq. (Q) with the dispersed frequency given by Eqs. 
d^-]^) are given in the Appendix. The final expression for the Bogoliubov coefficient |/3„p 
is: 

sinh^(27rfi_) +T{ko,B) 



\Pn\ 



sinh^(27rfi+) - sinh^(27rfi_) 



;ii) 



where Cli = Qi/n, Cl± = (fi""* ± fi*")/2, and the deviation function r{kQ,B) that contains 
the departure from thermality in the spectrum is 



,7r 



r(A;o, B) = cosh\-V4:Be~'=o _ i) . 



(12) 



When B = 0, for f2*" > fi""*, then it is clear from Eq. (11) that the spectrum of created 
particles is nearly thermal to high accuracy, 



Pn 



ttr 



-inVC 



(13) 



The function r(A;o, B) represents the deviation of the spectrum from thermal behavior due 
to the non-linearities at low momentum. Therefore, the amplitude of the power spectrum, 
Sffik), will be modified by r(fco, B) due to the non-linear dispersion function introduced at 
around Xq < 1. 



In de Sitter space, the Bogoliubov coefficients would not depend on k except their 
dependence in the bump parameters through the deviation function r(A;o) B). This function 
represents the departure from thermahty in the particle creation number, |/3nP and it 



confirms B. L. Hu idea [^ that near thermal radiance can be characterized by departure 
from exponential scaling. It is straightforward to derive the CMB power spectrum, P{n), 
analytically from (the exact solution for) the Bogoliubov coefficients an, Pn |]T^ 



^H = 7j^|-r^|/5n + «nP- (14) 

The deviation of the spectrum from scale-invariance in this class of models depends on the 
parameters of large-distance curvature modifications, namely: the scale of modified long 
wavelength modes, kQ^, and the deviation amplitude B. 

The expression for the Bogoliubov coefficient and Eq. (13) indicate that: for a well- 
defined initial vacuum state^, the spectrum is insensitive to the nonlinear dispersion relation 
in the transplanckian regime (modifications of short-distance physics). The unusual CMB 
spectrum plotted in the next Section with CMBFAST, demonstrates that modifications of 
the large scale curvature of the universe produce a tilt due to the departure from scale- 
invariance, and therefore conflict with the observed CMBR spectrum. In general the tilt is 
enhanced for modifications at superhorizon scales (fco ^ -f^o) because it is the low energy 
modes that dominate the spectrum in the Bogoliubov coefficient. Although departure from 
scale invar iance is smaller at the last scattering horizon scale, H^Sy the range of deviation 
parameters is constrained by the amplitude of the ffist peak. The deviation introduced to 
the spectral index, Ug from higher energy modes (wavelengths shorter than the last scat- 
tering horizon fco > Hls) becomes negligible because high energy modes do not contribute 
significantly to the spectrum. However, the shorter wavelengths would correspond to the in- 
termediate FRW regime rather than the large distance scales, a regime which is scrutinized 
by direct observation. 



3 CMB Spectrum 

Recent Boomerang and MAXIMA-1 CMB experiments |2l|, ^ have, to high accuracy, 
constrained the cosmological parameters, derived from the family of inflationary adiabatic 
models, to: total energy density Qtot = 0.90 ± 0.15 and spectral index Us = 0.99 ± 0.09 at 
a 95% confidence level [^. The current data favors a universe with dark energy density 
^A = 0.7. 

In this part, we explore the cosmological consequences of the alternative model that 
was given in Section 2 (Fig. 1). CMB is the most difficult test of precision cosmology 
that these models should pass. This model contains no dark energy, Q\ = 0, however it 
describes a universe which at large distances has a modified curvature from the metric of 
the FRW universe at intermediate scale. In Fig.2 we show the CMB power spectra obtained 
for different representative values of the deviation parameters ko and B in the dispersion 



^The field is in an initial Bunch-Davies vacuum. 
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function Eqs. (^0). The conventional parameters that go in the input of CMBFAST 
are: {fitot, ^b, ^c, ^a), which stand for total energy density, baryonic, cold dark matter and 
the cosmological constant energy density respectively; and Ug which is the scalar spectral 
index. We modified the power spectrum amplitude S^{k) in the POWERSFLAT subroutine 
of CMBFAST, in order to contain the deviation from the thermal spectrum (for the exact 
calculation reported in Section 2). The modified perturbation amplitude S'fj{k) is expressed 
in terms of 6%, ko, B, where 5^ is the unmodified amplitude of the scale-invariant power 
spectrum, ko corresponds to the location-scale where the curvature is modified, and B 
measures the amplitude of deviation in the curvature at scale ko- 

The values of the conventional parameters were taken to be (1,0.03, 0.97,0) for all the 
deviation plots (// — V), but the deviation parameters in the 4 plots below in Fig. 2 are 
in respective order: 

0, 5 = 0, ^A = 0.7) 

10-6 hMpc-\ B = 2) 

0.05 hMpc-\ B = 2) 

5 hMpc'\ B = 2) 

0.05 hMpc-\ B = 2.5) 
All plots were normalized to COBE. Shown for comparison is also plot / corresponding 
to the conventional CMB spectrum with ^a = 0.7. As it can be seen from the plots in Fig. 
2, there are distinct features of the CMB spectra corresponding to the dispersion function 
in comparison to the standard spectrum obtained for (A)CDM models. 

There is an overall tilt produced in the spectrum which signals departure from the 
scale invariance. This tilt is a function of the amplitude and scale of the modification, 
ko,B, introduced in Sect. 2 (Eq.|lT]), such that it increases for low values of the deviation 
momentum scale ko and large deviation amplitude B. Let us consider the 3 regimes into 
which the curvature modifications can be introduced: 

(1) Modifications at superHubble scales {ko < Ho). The departure from scale invariance 
is the strongest because the low energy modes dominate the spectrum, (II) in Fig. 2. 
Models predicting curvature modifications in regime (1) quite likely are ruled out due to a 
strongly tilted spectrum. 

(2) Modifications in the distance range between the current horizon Ho and last scat- 
tering horizon scale H^s {Ho < ko < His) ■ For this range, the tilt is less pronounced then 
in regime (1). The main constraint comes from the tilt and it tightly limits the amplitude 
of deviation in the first peak. For modifications around the last scattering horizon scale, 
ko ~ Hls the departure from scale invariance is vanishing, therefore the constraints are 
relaxed. However, even in this case the parameter B is tightly constrained to deviation by 
less than 10%, in order for the amplitude of the first peak Ai to be in the allowed range 
of 4500-5500 fxK'^ [0, ^ . In Fig. || we show the CMB spectra for these tuned values of 



ko, B; for comparison we also plotted the CMB for a value oi B = 2.5, which is outside the 
allowed range. 

(3) Modifications at distances shorter than the last scattering horizon {ko > His)- As 
we approach higher energy modes, the effect of the modification in the tilt of the spec- 
trum is suppressed, therefore the departure from the conventional spectrum is vanishing. 
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5000 




Figure 2: The CMB spectrum corresponding to our model for three different representative values of the 
deviation parameters kg and fixed deviation amplitude B = 2, II — IV. Shown is also the case of a larger 
amplitude B = 2.5 at scale fco « H^s (V). The CMB plots were obtain using CMB FAST and they were 
normalised to COBE. 

Nevertheless, these length scales do not correspond to large distances anymore, instead 
they are in the intermediate regime of FRW Universe. Thus the possibility of curvature 
modifications at such scales (galactic and intergalactic) is ruled out by direct observation 
up to very short distances (less than 1 mm). Clearly, there is no tilt or departure from the 
conventional CMB produced in the limit of modifications of very short distance physics, 
(very high momenta /cq ~^ oo? i-G-? transplanckian regime). 

The claim of the model was to "offer an alternative explanation" to the SNla data, 
namely: either the conventional Friedmann equation is modified or; the light of the SNla 
passing through regions of modified curvature would get deflected, and therefore when 
received by us would appear as if indicating an accelerating universe. Although this alter- 
native approach to the SNla data might be theoretically appealing, we conclude that the 
CMB data tightly constrains it and makes it unlikely to bear any resemblance to reality. 
The method used in this work, can also be adopted to check if higher dimensional models 
that predict modified gravity at large scales and modified equations for the perturbations 
□ [i i i i satisfy the CMB constraints. 



^These models naturally modify the curvature around horizon and Planck length-scales due to the 
higher dimensional gravity effects that switch on at very large or very short distances, but nevertheless 
with contributions from higher graviton excitations suppressed 124] . 



4 Summary 

In this work we investigated claims that a modified large- distance curvature may offer an 
alternative explanation for the SNla data. To check these claims, we studied the sensitivity 
of CMB spectrum to the whole range of modes, < fc < cxd, when short and large distance 
regimes are modified. 

In IP] it was noticed that a modified curvature of the universe at large distance (when 
traced back at the time of infiation^ gives rise to a dispersed frequency for the cosmic 
perturbations. The field is minimally coupled to the curvature thus its propagation feels 
the modifications in the background geometry. We adopted the method of Ref. ^ in order 
to find out the effects of curvature deviations on the current astrophysical observables. 

The role of a modified curvature of the universe at large distances on the infiationary 
metric perturbations was analytically described by a family of dispersion relations. The 
modification modulates the generalized frequencies of the infiationary perturbation modes 
at small values of the momenta k by departing from linearity around some certain small 
momenta k = k^ {ko < Mp) with a deviation amplitude B. The nonlinear feature of the 
dispersion relations, at small momenta k^ and in the transplanckian regime, tracks the 
curvature deviations at large and short distances, from the conventional FRW universe of 
intermediate scales. One of the parameters {k\ < k^,), in this class of dispersion functions 



was constrained in order to satisfy the Starobinsky bound for negligible backreaction [jT4|. 

The analytical expression for the CMBR spectrum (Sect. 2), as well as the CMBFAST 
plots of this class of models, deviate from the black body scale invariant spectrum. The 
deviation function r{kQ,B), given in Sect. 2 and the Appendix, which measures depar- 
ture from the scale-invariant spectrum (deviation from thermality in the the Bogoliubov 
coefficient), depends on two free parameters, the scale ko and the amplitude of the curva- 
ture modifications B. The tilt produced in the spectrum due to T{kQ, B) is present for all 
modification scales k^ < Mp (these values of the physical momenta correspond to the time 
of infiation). The tilt is less pronounced for scale modifications corresponding to length 
scales less or equal to the horizon of the last scattering surface, and in this case, the main 
constraint comes from the modifications to the amplitude of the first acoustic peak and the 
fact that curvature modifications in the intermediate FRW Universe scales are under direct 
observation. It remains interesting to answer why the only curvature modifications that for 
a small range of k^ and B can reconcile with the conventional CMB spectrum are allowed 
only around the last scattering H^s scales. 

The scale and amplitude of the deviations from the conventional spectrum, are severely 
constrained from the observed CMB spectrum to be within 10% of the scale and amplitude 
of the first peak. Although it is counterintuitive, since large distance would correspond to 
low energy theories, our results indicate that any modifications in the large scale curvature 
of the universe, is tightly constrained from CMB data, to a very small range of deviations 
from the curvature of the intermediate FRW universe. Perhaps, there is a natural way 
that would explain such a universe with an FRW spacetime at intermediate and very large 



*It should be noted that in the case of higher dimensional multigravity ^ ||, |^, |[ |2^, it is not clear 
how the metric perturbation equations are modified. 



distances but with small curvature deviations around His, without the need to appeal to 
fine-tuning. But if not, then theoretical cosmological models would have to account for the 
negative pressure dark energy of the universe. 

The analysis of the sensitivity of the CMB spectrum for the whole range of modes in 
a de Sitter background space, with modifications in the short and large distance physics, 
reveal: the spectrum is insensitive to the details of short- distance physics and the cuttoff 
scale kc (the transplanckian regime) only for an initial adiabatic vacuum state; the scale- 
invariance of the spectrum and the amplitude of the first acoustic peak are very sensitive to 
modifications of large distance physics (low momentum modes); the spectrum is also highly 
sensitive to the choice of the initial condition^ In our class of dispersion models, the 
initial vacuum state is well-defined since the background [QnivY) goes asymptotically fiat 



at early times (Bunch-Davis vacuum [^). The CMB spectrum for this class of models is 
indeed insensitive to short distance modifications, as it can be checked by taking the limit 
when large scale modification parameters ko, B go to zero, in which case the conventional 
scale-invariant spectrum is recovered. Therefore all the features observed in Fig. 0, are due 
to large-scale curvature modifications only. 
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5 Appendix 

The family of dispersion functions we used in Sect. 2 to model the deviation of the curvature 
at large and short distances is given by: 

F[k]^ = (e-kl)Vo{x,xo)+k^Vi{x-xo) + kl, (15) 

^i(---o) = -^ (i + e(-o))^ ' ('') 

where[3 x = k/kc, xq = ko/kc; kc = Mp is the cutoff scale, ko <C kc [xq <^ 1) is the 
value at which we deviate from linearity at low momentum, and the amplitude of the 
"bump /deviation" is controlled by B (see Fig. |l|). The parameter fci = ni/a{ri) gives the 
asymptotic constant value at initial time for the frequency in the limit {k ^ kc), i.e., in 
the transplanckian regime. On the other hand, in order to ensure the linear behavior at 



^It has been argued by many authors |^, |4| |T^ that the adiabatic vacuum is the right choice for the 
initial conditions. Even for the same dispersion model, a different choice for the initial vacuum state will 
clearly result in a different particle spectrum, therefore one has to be careful to distinguish if the features 
observed in the CMB spectrum are signatures of new physics or only of the choice of initial conditions. 
"'^"The momentum k has been shifted by kc such that fin ~ x for small positive values, x <t^ 1. 



very low values of the momenta, a; <^ 1, we impose the following constraints for any value 
of the deviation parameters Xq and B: 



Vo(x<l,Xo) ~ 1, 1/o(a; < l.xo)" ~ 0, 
Vl(x<l,a;o) ~ 0, 



(18) 
(19) 



where prime denotes derivative with respect to the physical momentum k. The generalised 
comoving frequency fln{v) is then given by: 



rinivf 



a{r])^F[n/a{ri)f 
C 



(n — n-, 



E e" 



— n 



B e 



(l + e(^-^o))2 



1 + e^ (l + e^)(l + e(^-^o)) 
with n = a{ri)k, a{rj) = —\ric\/ri during de Sitter inflation {\ric\ = l/H{ric)), and 

k 



+ nl, (20) 



n 



X 



-V- 



kc \Vc\kc 

The generalised frequency VLnirj) goes to constant values at r/ ^ ±oo, such that: 



fin(^) 

fin(^) 






^n=ni. 



^r = ^nl + C{n^ - nl) 



(21) 



(22) 
(23) 



with r2°"* ~ \fCn when n^ ni. 

Under the change of variables r] ^ u = exp{dnri), where d = l/{\ric\kc), the scalar wave 
equation (||) for the mode /i„ becomes: 



dl + -d^ + V{u) 
u 



/i„ = 



where: 



and, 



V{u) = t) + 



C 



+ 



E 



B 



u{l + u) m(m + 1)(7o + m) m(m + 7o) 



2 ' 



(24) 



(25) 



C = Clin" - nl)/{n df] , E = E[{n^ - nl)/{n df] 
D = {m/ndf, 
B = B/d^ , 



and 7o = exp{~-kQ / kc) ■ Eq. (|24D is exactly solvable in terms of the Riemann generalised 
hypergeometric functions [^ with the constraint E = C/{1 — 70), 



00 —70 
/in oc P I a c b u 
a* c* h* 



(26) 
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As explained in Section 2, because of the asymptotic behavior of Qniv)^ the initial 
and final vacua are well defined and the mode functions /i„ behave as plane waves in the 
asymptotic limits t] — > =1=00. The exact solution which matches this asymptotic behavior is 
then given by: 



;x™(r7) = N'^'iuyiu + 70)^ 2Fi[a + b + c,a + b + c*,l + a - a* 



where A^*" is a normalization constant, and 



u 



7o 



(27) 



a 
b 



-ifi*" = -i^D 



-(l + Vl-4i?e--o^ , 



-in"^^ = -id D + c 



(28) 
(29) 

(30) 



At late times the solution becomes a squeezed state by mixing of positive and negative 
frequencies: 



/ir(^) 



+ 



A^°"* (m)"(m + 7o)''x 

r(l + a-a*)r(c*-c) 
T{a + b + c*)r(l -a* -b-c) 
r(l+a-a*)r(c* -c) 



2Fi[a + b + c,a* + b + c,l + c-c*,-—] 

u 



T{a + b + c)r(l -a* -b-c* 



,Fi[a + b + c* ,a* + b + c* ,1 + c* - c, - 



70n 



u 



out 



OL. 



^ g-if^r*^ _^ 



/5„ 



,+^f^r'^ 



2fi°"* 



2fi°"* 



(31) 



with |/9riP being the Bogoliubov coefficient equal to the particle creation number per mode 
n and Vti = Qi/n. Using the linear transformation properties of hypergeometric functions 



p6| , we find that 



Pn 



Olr, 



sinh^(27rfi_) +r(/to,5) 
sinh^(27rfi+) + r(fco,B) ' 



where. 



VLa 



Qout _^ Qin 



n^i 



nb 



(32) 



(33) 



and the deviation function r(A;o, B) is 



V{kQ,B) = cosh^(-V45e-^o - 1) 



(34) 



When 5 = and fi*" > fi""*, then it is clear from Eq. (^) that the spectrum of created 
particles is nearly thermal to high accuracy, 



Pn 



ttr 



-4nVC 



(35) 
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as expected in de Sitter expansion. However, when i? 7^ 0, 70 7^ 0, at the mode crossing 
time n = Ha{ri), we can write: 



an 



^-i-n-VC 



1 + 



T{ko,B) 
sinh^ 2nQ- 



1 + 



r{ko,B) 

sinh^ 2ttQ4 



(36) 



The expression in the squared bracket in the above equation contains the deviation from 
scale invariance. The deviation T{ko, B) is larger at low values of the momentum modifica- 
tion scale, xq ^ 1. On the other hand, V{kQ,B) is suppressed around large scales, xq ~ 1. 
The same results about the scale dependence of the deviation function were obtained by 
using CMBFAST code (Figs. 0). 
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